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Abstract. We develop a general framework for finding error estimates for 
' convection-diffusion equations with nonlocal, nonlinear, and possibly degener- 

ate diffusion terms. The equations are nonlocal because they involve fractional 
diffusion operators that are generators of pure jump Levy processes (e.g. the 
fractional Laplacian). As an application, we derive continuous dependence es- 
timates on the nonlinearities and on the Levy measure of the diffusion term. 
Estimates of the rates of convergence for general nonlinear nonlocal vanishing 
viscosity approximations of scalar conservation laws then follow as a corollary. 
, Our results both cover, and extend to new equations, a large part of the known 

' error estimates in the literature. 

^' 

1. Introduction 
This paper is concerned with the following Cauchy problem: 

^ _ idMx,t)+dwifiu)){x,t)=a^[A{u{;tMx) in Qt :=K'^ X (0,T), 

00' \u{x,0) = uo{x), in R'^, 

00 ' 

' where u is the scalar unknown function, div denotes the divergence with respect to 

^SJ , (w.r.t.) X, and the operator £^ is defined for all (j) G C^{R'^) by 

O! (1.2) C''mx):= [ 4>ix + z)-(f>ix)-z-D^{x)l^,^<,dfi{z), 

where D(j) denotes the gradient of (j) w.r.t. x and l|z|<i — 1 for ^ 1 ^-nd = 
otherwise. Throughout the paper, the data (/, A, uo,/x) is assumed to satisfy the 
following assumptions: 

(1.3) / = (A, ...Jd)e W^i-°°(R,R^) with /(O) - 0, 

(1.4) A e W^^'^iR) is nondecreasing with A{0) = 0, 

(1.5) uo£ L'^{R'^)nL\R'^)nBV{m.'^), 
and 

(1.6) ^ is a. nonnegative Radon measure on R'^ \ {0} satisfying 



I 

JR" 



|zp A 1 d^(z) < oo, 

^\{0} 

where we use the notation a Ab = min{a, b}. The measure is a Levy measure. 



Remark 1.1. 

(1) Subtracting constants to / and A if necessary, there is no loss of generality 
in assuming that /(O) = and A{0) — 0. 
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(2) Our results also hold for locally Lipschitz-continuous nonlinearities / and A 
since solutions will be bounded; see Remark 2.3 for more details. 

(3) Assumption (1.6) and Taylor expansion reveals that jC^[4i\ is well-defined 
for e.g. bounded functions 0; 

\/:^mx)\<m^x\D^ct>{x + z)\ [ hzfd^,{z) + 2\\c|>\\L^ [ d^^) 

kl<l yO<|^|<l ^ J\z\>l 

where D'^(j) is the Hessian of </>. If in addition D'^<p is bounded on M.'^, then 
so is [(/)] . 

Under (1.6), £^ is the generator of a pure jump Levy process, and reversely, 
any pure jump Levy process has a generator of like (see e.g. [6, 54]). This 
class of diffusion processes contains e.g. the a-stable process whose generator is the 
fractional Laplacian - (- A)"^ with a € (0, 2). It can be defined for all (j) € C^{R'') 
via the Fourier transform as 

(-A)^ = .F-i (I • r.F(/.) , 

or in the form (1.2) with the following Levy measure (see e.g. [6, 32]): 

dz 

(1.7) dfx{z) = (up to a positive multiplicative constant). 

\z\ " 

Many other Levy processes/operators of practical interest can be found in e.g. 
[6, 24]. Under assumption (1.4), £''[A(-)] is an example of a nonlinear nonlocal 
diffusion operator. For recent studies of this and similar type of operators, we refer 
the reader to J8, 9, 15, 19, 27] and the references therein. 

Equation (1.1) appears in many different contexts such as overdriven gas detona- 
tions J22], mathematical finance [24], flow in porous media [27], radiation hydrody- 
namics [51, 52], and anomalous diffusion in semiconductor growth [57]. Equations 
of the form (1.1) constitute a large class of nonlinear degenerate parabolic integro- 
differential equations (integro-PDEs). Let us give some representative examples. 
When ^ = or ^ = 0, (1.1) is the well-known scalar conservation law (see e.g. [25] 
and references therein): 

(1.8) dtU + diYf{u) =0. 

When A{u) = u, (1.1) is the so-called Levy/fractal/fractional conservation law: 

(1.9) dtu + diY f{u)=jC^[u]. 

Equation (1.9) has been extensively studied since the nineties [1, 2, 3, 4, 5, 7, 10, 
11, 12, 16, 17, 20, 21, 26, 28, 29, 30, 31, 32, 35, 37, 38, 39, 40, 41, 44, 48, 49, 50]. 
When A is nonlinear, (1.1) can be seen as a generalization of the following classical 
convection-diffusion equation (possibly degenerate): 

(1.10) dtu + div/(u) = AA{u); 

see e.g. [13, 14, 18, 23] for precise references on (1.10). The case of nonlinear and 
nonlocal diffusions has been studied in [27] in the setting of nonlocal porous me- 
dia equations, and in [19] where a general L^-theory for (1.1) is developed along 
with connections to Hamilton-Jacobi-Bellman equations of stochastic control the- 
ory. Other interesting examples concern the class of nonsingiilar Levy measures 
satisfying /r<i^|o} '^l^i^) ^ +oo. In that case, £^ is a convolution operator and (1.1) 
can be seen as a generalization of Rosenau's models [42, 43, 47, 48, 55, 56] and 
nonlinear radiation hydrodynamics models [51] of the respective forms 

(1.11) dtU + div/(u) = * u — u, 

(1.12) dtu + dWf{u)=g^*A{u)-A{u), 



CONTINUOUS DEPENDENCE ESTIMATES FOR INTEGRO-PDES 



3 



where * denotes the convolution product w.r.t. x and G is nonnegative 

with fj^ag^{z)dz= 1. 

Most of the results on such nonlocal convection-diffusion equations concern Equa- 
tion (1.9) whose diffusion is linear. It is known that shocks can occur in finite 
time [4, 28, 42, 43, 44, 48, 55], that weak solutions can be nonunique [2], and that 
the Cauchy problem is well-posed with the notion of entropy solutions in the sense 
of Kruzhkov [1, 41, 47, 55]. Results on nonlinear nonlocal diffusions can be found 
in [51] where entropy solutions of (1.12) are studied. Very recently, the entropy 
solution theory has been extended in [19] to cover the full problem (1.1) for general 
singular Levy measures and nonlinear A. 

The purpose of the present paper is to develop an abstract framework for find- 
ing error estimates for entropy solutions of (1.1). As applications, we focus in this 
paper on continuous dependence estimates and convergence rates for vanishing vis- 
cosity approximations. We refer the reader to [13, 18, 23, 46] and the references 
therein for similar analysis on (1.10) and related local equations. As far as non- 
local equations are concerned, contimioiis dependence estimates for fully nonlinear 
integro-PDEs have already been derived in [36] in the context of viscosity solutions 
of Bellman- Isaacs equations; see also [32, 34. 36] for error estimates on nonlocal van- 
ishing viscosity approximations. To the best of our knowledge, there arc only a few 
results for nonlocal conservation laws. All the results we have found concern Equa- 
tions (1.9), (1.11) and (1.12) for which wc refer the reader to [1, 29, 32, 41, 47, 55]. A 
large part of these error estimates concern convergence rates for vanishing viscosity 
approximations. The only result we have found on continuous dependence estimates 
appears in [41]; it concerns Equation (1.9) in the case of self-adjoint Levy operators. 
To finish with the bibliography, let us also refer the reader to [20, 21, 26, 30, 51] 
for the related topic of error estimates for numerical approximations. 

Our main result is stated in Lemma 3.1, and it compares the entropy solution u 
of (1.1) with a general function v. Our main application consists in comparing u 
with the entropy solution v of 

^^_^3^ (dtV + divgiv)^£-^[B{v)], 

\v{x,0) = vo, 

where the data set {g, B, vq, u) is assumed to satisfy (1.3)- (1.6). We obtain explicit 
continuous dependence estimates on the data stated in Theorems 3.3 3.4. Let 
us recall that when i? = or = 0, (1.13) is the pure scalar conservation law 
in (1.8). Equation (1.1) can thus be seen as a nonlinear nonlocal vanishing viscosity 
approximation of (1.8) if A or /i vanishes. The rate of convergence is then obtained 
as a consequence of Theorems 3.3-3.4, see Theorem 3.7. 

It is natural to compare Theorems 3.3-3.4 and Theorem 3.7 with the known 
error estimates for Equations (1.9), (1.11) and (1.12). One can see that a quite 
important part of them are particular cases of our general results. We discuss 
this point in Section 3 by giving precise examples. Let us mention that we also 
give a simple example of Hamilton- Jacobi equations suggesting that Theorems 3.3- 
3.4 are in some sense the "conservation laws' versions" of the results in [36]; see 
Example 3.2. 

To finish, let us mention that in the case of fractional Laplacians of order a > 1, 
Theorems 3.3-3.4 can be improved by taking advantage of the homogeneity of the 
measures in (1.7). In order not to make this paper too long, this special case 
(including a < 1) is investigated in a second paper [3]. 

The rest of this paper is organized as follows. In Section 2 we list the notation 
used throughout the paper; we also recall the notion of entropy solution to (1.1). 
In Section 3, we state and discuss our main results. Sections 4-5 are devoted to the 
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proofs of our main results; Section 4 states some preliminary results on the nonlocal 
operator. For the reader's convenience, some more or less known technical results 
have been gathered in the Appendices. 

2. Preliminaries 

In this section we explain most of the notation used in the paper, and we give 
the definition of entropy solutions of (1.1) along with a well-posedness result. For 
the definitions of measures and BV-spaces, we refer to the books [33, 53] . 

2.1. Notation. 

2.1.1. Vectors, sets and functions. Throughout the paper d € N is a fixed dimen- 
sion, T > a fixed time, and {x,t) = {xi, . . . ,Xd,t) £ Qt ■= R'' x (0,T) is the 
generic space-time variable. For all a,b G R we let a A b := min{a, b}, a V 6 := 
max{a. b}. a+ := a V 0, and a~ := (—a) V 0. For all m G N. wc let • and | • | denote 
the Euclidean inner product and norm of M™, while for a matrix A G R™^™, we use 
the norm \A\ = max{Aw : w G M™, |u;| < 1}. We let -E := {-w gW^ : w G E}, 
and denote the characteristic function of the set E' by 1^;. 

By and we denote the spaces of infinitely differentiable functions and 
infinitely diff'erentiable functions with compact support. Moreover, for p G [1, +00], 
LP, W^'P, Lj^^, and V denote the Lebesgue and Sobolcv spaces, the locally inte- 
grable functions, and the Schwartz distributions respectively. Sometime we indicate 
range and domain of the functions, e.g. C^ (K'^M'*). 

The support oi u G V is denoted by suppu. The restriction of m to a set U is 
denoted by u^^j. By u * v we mean the convolution of two functions u = u{x,t) 
and V = v{x,t) w.r.t. to the space variable x. We let dtu, D^u, and D'^u denote 
the partial derivative in time, the; spatial gradient, and the spatial Hessian matrix 
of u respectively. If there is no confusion, we write D instead of D^- The derivative 
of a one variable function u is written u'. The same notation is also used for 
distributional derivatives. 

2.1.2. Radon measures. Let n he a nonnegative Radon measure on M'' \ {0}, i.e. a 
measure n : BRd\{o} ~^ [Oi +00] which is finite on compact sets and where BRd\{o} 
is the Borel cr-algebra of M."^ \ {0}. As usual, fi is extended to a complete measure 
to the (smallest) /x-completion of Bg^d^^Qy: 

Sm<*\{0} := {ECR'i\ {0} s.t. there is B, G Br.\{o} (* = 1,2) 
with B1QEQB2 and ^1(^2 \ Bi) = 0}. 

We say that u : K'' \ {0} ^ R is /i-mcasurablc (/i-integrable) if for each real inter- 
val I, u~^{I) G -Bgd^^Qj (if in addition it is integrable w.r.t. ji). For the Lebesgue 
measure we simply use the terminologies measurable, integrable, almost everywhere 
(a.c), etc. 

Throughout the paper the Lebesgue measure of R™ is denoted by Aw if w denotes 
the generic variable of R™. Its tensor product with is denoted by d/i(2) Aw; note 

that this is a well-defined nonnegative Radon measure on j^^]^'i\^{o})xiR^' since 
is cr-finite. 

Given another nonnegative Radon measure v on \ {0}, the total variation 
oi ^ — V is denoted by |/i — z/l. The positive and negative parts are 
denoted by (/x — z/)^ respectively. 
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2.1.3. BV-spaces. The space BV{M.''') of functions with bounded variation on R'^ 
is defined as the space of m e Ll^^{R'^) such that Du is a (finite variation Radon) 
measure Du : Bgd R'^. Let us recall that its total variation \Du\ satisfies for 
all B G B^d, 

(2.1) 

\Du\{B)= inf supl / udiv^dx: ^gC?°{R'^, R'^),\(f)\<l, supple u\ . 

BCU open [^Jjjd J 
The i?y-scmi-norni of u is defined as |?i|BV'(K<i) := |-Du|(M'') < +oo . 

2.1.4. Functions in L°"nC{L^)nL°"{BV). Given u e L~(Qt) n C([0, T]; L^) we 
denote by u{t) the function u{-,t) if there is no confusion. The C([0, T]; I/^)-norm 
of u is 

l|w||c([o,T];Li) := max |||u(i)||ii(Kd) : t € [0,T]|. 
The modulus of continuity (in time) of u € C([0, T]; L^) is denoted by 

a;„((5) :=max{l|M(t) - w(.s)|lii(R<i) : t,s€ [0,T] s.t. \t - s\ < 6} , 6 > 0. 
Throughout the paper wc say that 

u e (Qt) n C([0, T] ; L^) n L°° (0, T; BV) 
if for ah t e [0, T], zi(<) G ^^(M''), and if the i°°(0, T; BF)-scmi-norm of u, 

\u\l'^{o,T;BV) ■■= sup ||u(t)|i3v^(R<i) : t e [0,T]| < 00. 
The L^{0,T; By)-semi-norm of u is defined as 

|'"Ui(o,T;By) := / |'"(*)|By(R<i) di. 
Jo 

Note that it is standard that t \u{t)\Bv(R'i) is measurable since it is lower semi- 
continuous; see e.g. [33]. 

2.2. Entropy formulation and well-posedness. Let us recall the formal com- 
putations leading to the entropy formulation of (1.1). First we split into 3 

parts: 

(2.2) C'^mx)^C>^mx)+diy{b>^<P){x)+C^'^mx) 
for (p e C^{R'^), r > 0, and a; e R'', where 

(2.3) C'^[ct>]{x):= [ (l>{x + z)-cl>{x)-z-Dct>{x)l\,i<,dii{z), 

Jo<\z\<r 

(2.4) K~-[ ^l|z|<idM(^), 

J\z\>r 

(2.5) £^''-[(^](a;) := / <P{x + z) - <f>{x) dtx{z) . 

J\z\>r 

Consider then the Kruzhkov [45] entropies | • — A;|, fc € K, and entropy fluxes 

(2.6) qf{u, k) := sgn {u - k) {f{u) - f{k)) e 

where we always use the following everywhere representative of the sign function: 

2.7) sgn(M :=<^ 

10 if u = 0. 

By (1.4) it is readily seen that for all u,k & R, 

(2.8) sgn (u - k) {A{u) - A{k)) = \A{u) - A{k)\, 
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and we formally deduce from (2.2), (2.8), and the nonnegativity of ji that 
sgn{u-k)C''[A{u)] 

< £^[|A(u) - A{k)\] + div {V^ \A{u) - A{k)\) + sgn {u - k) £''''-[A(«)]. 

Let u be a solution of (1.1), and multiply (1.1) by sgn (it — k). Formal compu- 
tations then reveals that 

dt\u -k\+ div (g/(u, k) - b'f \A{u) - ^(fc)|) 

< C>^[\A{u) - A{k)\] + sgn {u - k) C'^^'^lAiu)]. 

The entropy formulation in Definition 2.1 below consists in asking that u satisfies 
this inequality for all entropy-flux pairs (i.e. for all fc e ffi) and all r > 0. Roughly 

speaking one can give a sense to sgn [u — k) £^^''[A(u)] for bounded discontinuous u 
thanks to (1.6). But since ji may be singular at z = 0, see Remark 1.1 (3), the 
other terms have to be interpreted in the sense of distributions: Multiply by test 
functions (f) and integrate by parts to move singular operators onto test functions. 
For the nonlocal terms this can be done by change of variables: First take {z, x, t) — >■ 
{—z, X, t) to see (formally) that 

/ (j>dYv {hi^\A{u)- A{k)\)Axdt= [ D(j> ■ h^i:' \A{u) - A{k)\dxdt, 

JQt JQt 

where /x* is the Levy measure (i.e. it satisfles (1.6)) defined by 
(2.9) ii*{B):=n{-B) for all B€B^d\{o^. 

In view of (2.3), we can take {z, x, t) {—z, x + z,t) to find that 

[ (l)Ci^[\A{u)- A{k)\]dxdt= [ \A{u) - A{k)\jCf[(l>]dxdt. 

JQt JQt 

This leads to the following definition introduced in [19]. 

Definition 2.1. (Entropy solutions) Assume (1.3)-(1.6). We say that a func- 
tion u e L°°{Qt) n C ([0,T];L^) is an entropy solution of (1.1) provided that for 
all A; e R, allr > 0, and all nonnegative (j) e C^(K''+^), 



(2.10) 



/ \u-k\ dt<l) + (qf{u, k) + bf \A{u) - A{k)\) ■ Dcjidxdt 

JQt ^ ' 

-h / |A(M)-A(fc)|rf [(?i]-Fsgn(u-A;)£'''''[74(u)]0da;dt 

JQt 

- \u{x,T)-k\(l){x,T)dx+ \uo{x) - k\(l){x,0)dx>0. 
Jm.'' Jm.'' 

Remark 2.1. 

(1) Under assumptions (1.3)-(1.6), the entropy inequality (2.10) is well-defined 
independently of the a.e. representative of u. To see this note that n* ob- 
viously satisfies (1.6), and hence it easily follows that £{f [(j)] € {M.'^'^^) . 
Since sgn(w-fc), qf {u,k), and A{u) belong to L°° by (2.7) and (1.3)-(1.4), 
it is then clear that all terms in (2.10) are well-defined except possibly the 
£''^''-tcrm. Here it may look like we are integrating Lebcsguc measurable 
functions w.r.t. a Radon measure /x. However, the integrand does have 
the right measur ability /integr ability by Lemma 4.2. We therefore find that 
since A{u) belongs to C{[0,T]; L^), so does also C'^'^[A{u)] and we arc done. 

(2) In the definition of entropy solutions, it is possible to consider functions u 
only defined for a.e. t G [0, T] by taking test functions with compact support 
in Qt and adding an explicit initial condition, see e.g. [19]. 
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(3) One can check that classical solutions are entropy solutions, thus justifying 
the formal computations leading to Definition 2.1. Moreover entropy so- 
lution are weak solutions and hence smooth entropy solutions are classical 
solutions. We refer the reader to [19] for the proofs. 

Here is a well-posedness result from [19]. 

Theorem 2.2. (Well-posedness) Assume (1.3)-(1.6). There exists a unique en- 
tropy solution u of (1.1). This entropy solution belongs to L°°{QT)r\C ([0, T]; L^) n 
L°° {0,T;BV) and 

{||w||l~(Qt) < hoWL^-iR"), 
ll"llc([0,T];Li) < II"o||l1(R''), 
\u\l°°{0,T;BV) < |uo|BV(R'i)- 

Moreover, if v is the entropy solution of (1.1) with v{0) = vo for another initial 
data Vo satisfying (1.5), then 

(2-12) \\U - v|lc([0,T];Li) < jjwo - ^0 ||z,i(M<i) • 

Remark 2.3. By the L°°-estimate in (2.11), all the results of this paper also holds 
for locally Lipschitz-continuous nonlinearities (/, A). Simply replace the data (/, A) 
by if, A) 1[_M,M] with M := ||uo||z,oc(M<i). 



3. Main results 

Our first main result is a Kuznetsov type of lemma that measures the distance 
between the entropy solution u of (1.1) and an arbitrary function v. 
Let €,6>0 and (j)"'^ e C°°{Qt) be the test function 

(3.1) (l)''\x, t, y, s) := Os{t - s) 0,{x - y), 

where 9i{t) := and 9f:{x) := ^^d(f) are, respectively, time and space 

approximate units with kernel On with n = 1 and n = d satisfying 

(3.2) 9neC^{R"), en>0, supp4 C {|a;| < 1}, and / ^„(x)da;=l. 

Also recall that u;„((5) is the modulus of continuity of u e C ([0,T];X^). 

Lemma 3.1 (Kuznetsov type Lemma). Assume (1.3)-(1.6). Let u be the entropy 
solution of (1.1) andv e L°°{QT)nC ([0, T]; L^) with v{0) = vq. Then for allr > 0, 
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e > 0, andO<S <T, 



(3.3) 

\\u{T)-v{T)\\m^,) 



< Iho - wo||z,i(R<J) +^Cg\uo\BV{M<^) +2a;„(5) Vw^(^) 




ff ^ (qfiv{x, t),u{y, s)) + bf \A{v{x, t)) - A{u{y, s))|) • D^(t>''\x, t, y, s) dw 



+ [[ \A{v{x,t)) - A{u{y,s))\jrf[<l>''\x,t,-,s)]{y)dw 



JjQl 



[[ sgn {v{x, t) - uiy, s)) £''''-[^(w(-, s))](y) <P''^{x, t, y, s) dw 



+ [[ \v{x,T) - u{y,s)\^''\x,T,y,s)dxdyds 




R-'xQt 



\vo{x) - u{y,s)\^'''^{x,0,y,s)dxdyds 



where dw := dxdtdyds, and := 2 \x\6d{x) dx. 



Remark 3.2. 

(1) The error in time only depends on the moduU of continuity of u and v 
at t = and t = T. Here we simply take the global-in-time moduli of 
continuity uju{S) and ujy{S), since this is sufficient in our settings. 

(2) When A = or /x = this lemma reduces to the well-known Kuznetsov 
lemma [46] for multidimensional scalar conservation laws. 

(3) Notice that the £{f -term vanishes when r ^ 0, sec Lemma 4.6. 

(4) Lemma 3.1 has many applications. In this paper and in [3] we focus on con- 
tinuous dependence results and error estimates for the vanishing viscosity 
method. In a future paper, we will use the lemma to obtain error estimates 
for numerical approximations of (1.1). 

In this paper we apply Lemma 3.1 to compare the entropy solution u of (1.1) 
with the entropy solution v of (1.13). This is our second main result , and we 

present it in the two theorems below. The first focuses on the dependence on the 
nonlincaritics (with /i = v) and the second one on the Levy measure (with A = B). 

Theorem 3.3. (Continuous dependence on the nonlinearities) Let u and v be the 
entropy solutions of (1.1) and (1.13) respectively with data sets {f, A,uo, n) and 
{g, B, vo, V = ^J') satisfying (1.3)-(1.6). Then for all r2 > ri > 0, 



U - v\\c([0,T];L'^) < \\U0 - VoW L'^ (Rd) + |wo|sy(M<i) T \\f' - g'\\L' 



(3.4) 
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where ca = |^ • 

Theorem 3.4. (Continuous dependence on the Levy measure) Let u and v be 
the entropy solutions of (1.1) and (1.13) respectively with data sets {f,A,uo,iJ.) 
and {g, B = A, vq, v) satisfying (1.3)-(1.6). Then for all r2> ri> 0, 

11^ - 'i^llc([o,T];Z,i) < ||wo - i^oIIl^K'') + ko|sy(Md) T II/' - g'WL"' 



(3.5) 



+ |Mo|By(R^) Jcdrp'lli^(K) / \z\'^ - v\{z) 

y ^0<lzl<n 

+ |wo|By(Md) 2^ m'||L~(R) / 

J n < 



/0<|z|<ri 

\z\d\^i-v\{z) 



L<|zj<r2 



+ I"o|bv'(r<') T II^'I|l~(r) 



/ zd{n-u){z) 

^riAl<|z|<riVl 



■T\\A'\\i^^(K.) / \\uo{-+z)-UQ\\Li(m.d)d\ii-u\{z), 

J\z\>r2 



where Cd 



d+l ■ 



From Theorems 3.3 and 3.4 we can easily find a general continuous dependence 
estimate when both A and /x are difi'erent from B and i^, respectively. E.g. we can 
take an intermediate solution w oi Wt + div f{w) = £'^[B{w)] and w{0) = uq, and 
use the triangle inequality. Using this idea we can show that the following estimates 

always have to hold: 

Corollary 3.5. Let u and v he the entropy solutions of (1.1) and (1.13) respectively 
with data sets {f, A,U(), fi) and {g,B,vo,u) satisfying (1.3)-(1.6). Then 

\\u - f llc([O.T];Li) < I|W0 - 'fo||Li(R'') + Wo\ BV (R^) T\\f' - 

^^■^^ + C (T^ V T) (^^\\A' - i?'|U=>o(M) + ^1^^^^^^ I^P A 1 d|/x - u\{z)^ 

where C only depends on d and the data. Moreover, if in addition 



\z\ A 1 dii{z) + \z\Al du{z) < oo, 

R<i\{0} yR<*\{0} 



then we have the better estimate 

\U - v\\c{lO,T];L^) < \\U0 " Vo\\Li(Rd^ + \Uo\bv{R'') T \\f' - ^'IIl" 

(3.7) 

i- I — ||L~(R) -I- 

/R<i\{0} 

where C only depends on d and the data. 



+ CT (\\A' - B'IUoo(M) + / 1^1 A 1 d|M - v\{z)\ , 

\ Jr-'\{o} / 



Outline of proof . To prove (3.6), we use Theorems 3.3 and 3.4 with n = 1 = r2 

and the triangle inequality. Wc also use estimates like |a — fej < y^jal + |6|-\/|a — b\, 
~ '^l ^ ImI + I'^l etc. To prove (3.7) we take n = and r2 = 1. □ 

Remark 3.6. 

(1) All these estimates hold for arbitrary Levy measures /z, i/ and even for 
strongly degenerate diffusions where A, B may vanish on large sets. They 
are consistent (at least for the — i^j term) with general results for nonlocal 
Hamilton- Jacobi-Bellman equations in [36] . When /x, v have the special 
form (1.7) (with possibly different a's), then it is possible to use the extra 
symmetry and homogeneity properties to obtain better estimates, see [3]. 
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(3. 



(2) The optimal choice of the ri, r2 depends on the behavior of the Levy mea- 
sures at zero and infinity, see the discussion above and at the end of this 
section for more details. 

(3) In the special case of symmetric Levy measures (/x = /x*), the terms corre- 
sponding to the cutting nAl < \z\ <riVl disappear. 

Let us now consider the nonlocal vanishing viscosity problem 

'dtu" + div/(«^) = eCi^[A{u<')], 
M«(0) = uo, 

i.e. problem (1.8) with a perturbation term e C^[A{u'^y\. When e > tend to zero, 
is expected to converge toward the solution u of (1.8). As an immediate application 
of Theorem 3.3 or 3.4, we have the following result: 

Theorem 3.7 (Vanishing viscosity). Assume (1.3)-(1.6). Let u and he the 

entropy solutions of (1.8) and (3.8) respectively. Then 

w1lc([o,Tl;Li) < C" min 

r2>ri>0 

(3.9) 

+ Te / \z\An{z)+ I zdniz) + I d/i(2;) 

.Jri<\z\<r2 ^riAl<|2|<riVl J \z\>r2 

where C only depends on d, \\2iQ 

Outline of proof. Note that u can be seen as the entropy solution of (1.1) with A = 
and n as Levy measure. Hence we can estimate w'^||c([o,t];L1) from Theorem 3.3. 
The error coming from the difFcrcncc of the derivatives of the nonlinearities is equal 
to e IIA'll^oo(jj'). Inequality (3.9) then follows from (3.4). □ 

Corollary 3.8. Assume (1.3)-(1.6). Let u andu"^ be the entropy solutions of (1.8) 
and (3.8) respectively. Then 

\\u-u'\\ci[o,T];L^)<CiT^WT)eK 
where C only depends on d and the data. Moreover, if in addition 

\z\ A 1 diJ,{z) < 00, 




/ 



then we have the better estimate 

11^*- "1lc([0,T];Li) < CTe, 
where C depends on d and the data. 

This corollary follows immediately from Theorem 3.7 or Corollary 3.5. 

Remark 3.9. 

(1) Our estimates are just as good or better than the standard 0{ei) estimate 
for the classical vanishing viscosity method ((1.10) with A{u) = en). 

(2) Our estimates hold for arbitrary Levy measures /x and even for strongly 
degenerate diffusions where A may vanish on a large set! This is consistent 
with general results for nonlocal Hamilton- Jacobi-BcUman equations [36]. 

(3) If the solutions are smoother, it is possible to obtain better estimates. Also 
if /It is as in (1.7), the additional symmetry and homogeneity can be used to 
obtain better estimates which can be proven to be optimal. See Example 
3.3 below. 

(4) Corollary 3.8 contains less information than Theorem 3.7 and is not strong 
enough to get optimal results in all cases, e.g. in Example 3.3 with a > 1. 
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(5) The error estimates above trivially also holds for the more general vanishing 
viscosity equation 

jdtu'' + div/Cu-^) = C''[B{u'')] +eCi^[A{u% 
\u'{0) = uo. 

Further discussion. We now make a more precise comparison of the results above 
with known estimates from the literature. Wc begin with continuous dependence 
estimates and finish with convergence rates for vanishing viscosity approximations. 

Let u and v denote the entropy solutions of (1.1) and (1.13), respectively. To 
simplify, we take the same data sets (/, j4,?/o) = {fj,B,i;a) and we only allow the 
Levy measures n and v to be different. We also let C denote a constant only 
depending on T, d and the data. 

Example 3.1. Let us consider Equation (1.9), i.e. A{u) = u. Let us also consider 
the class of Levy operators satisfying 

For such kind of equations, the following continuous dependence estimate on the 
Levy measure has been established in [41]: 

\\u-v\\c(xo,nLr)<C [ \zfd\fx-u\{z) + c[ \z\d\i,-t^\{z). 

^0<|z|<l ^|2|>1 

This estimate follows from Theorem 3.4 by taking n = 1 and r2 = +oo in (3.5). 

Example 3.2. Consider the following one-dimensional Hamilton-Jacobi-Bellman 
equation 

Ut + fiU,)=C^[U] 

with initial data Uo{x) := f^^uo{y) dy. This particular equation is related to 
the nonlocal conservation law (1.8), its solution U{x,t) = J^^u{y,t) dy where u 
solves (1.8), see [19]. It is also an example of an integro-PDE for which the general 
theory of [36] applies, and this theory allows us to establish the following continuous 
dependence estimate on the Levy measure: 



sup \U-V\<CJ Izl"^ Ald\ii-iy\{z), 

Mx[0,T] ]] Jr\{0} 

where V{x,t) := v{y,t)dy. (This result is a version of Theorem 4.1 (in [36]) 
which follows from Theorem 3.1 by setting Po, ■ ■ ■ ,Pa,Ps = and p = |^;[ A 1 in 
(AO)). Since 

sup < ||u-t;[|c([o,T];Li), 

Rx[0,T] 

this estimate also follows from (3.6) in Corollary 3.5 when (A, /, uq) = iB,g,VQ). 

Let us now compare Theorem 3.7 with known convergence rates. We keep the 
same notation for u, and 0{-) as in Theorem 3.7. 

Example 3.3. Let us consider the case where A{u'^) = and = — (— A)f , 
a e (0,2). Then the following optimal rates have been derived in [1, 29]: 



(3-10) \\U - u''\\cilO,T];L^) = < 



0[f--j ifa>l, 
0(e[lne[) if a = 1, 
0{e) ifa<l. 
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Let US explain how these results can be deduced from (3.9). First we use (1.7) to 
explicitly compute the integrals in (3.9) and obtain 

\\U - 'U^||(7l'[n.Tl:7-_,11 = O 




C([0,T];Li) 

We then deduce (3.10) by taking ri = and r2 = +oo if a > 1, ri = e and r2 — I 
if a = 1, and ri = and r2 = 1 if a < 1. 

Example 3.4. Let us consider the class of Levy operators where d^{z) = gix{z) dz 
for < .9^ G L^(R'^) such that ^^^gij_[z)dz = 1. This corresponds to problem 
(1.11)-(1.12) since wc may write 

[u'^] = g^* u'^ — (* is the convolution in space) . 

The following optimal rate of convergence has been derived in [47, 55] : 

\\u-u''\\ci[0,T];L^) =0{€). 

This result also follows from (3.9) by taking n = r2 = 0. 

4. Auxiliary results concerning £^ 

Before proving our main results in the next section, we state and prove some 
results concerning that will be needed. 



Lemma 4.1. Assume (1.6) and r > 0. Then for all (p & C[ 

WrW\\\LHR'^)< / kPd/^(^) ll0lllV2.i(R<i)- 

JO<\z\<r 

Proof. Recall that fj. is cr-finite as nonnegative Radon measure, hence the product 
measure of n and the Lebesgue measure is a well-defined nonnegative Radon mea- 
sure for which Fubini's theorem applies. Using in addition Taylor's formula with 
integral remainder, we find that 

\\j!:m\LHm'^) 

dx 



= / (j){x + z) - (j){x) - z ■ D(l){x)li^i<idfi{z) 

Jr<' Jo<\z\<r 

< [ [ 1^1^(1 -r) / \D^<l){x + Tz)\dxdTdn{z) 

Jo<\z\<rAlJo J'R'i 

+ [ I \A I \D<t){x+Tz)\dxdTdlji{z), 

Jl<\z\<rylJ<i Jr<^ 

< \z\'^dn{z)\\D^(f)\\Li^T^d^^4x4) 
Jo<\z\<rAl 

+ / \z\ dn{z) ||-D^||Li(R<i,K<i) (by the change of variable x x + t z), 

< / \zfdi^{z)\\<p\\w^.Hm- 

JO<\z\<r 



l0<\z\<r 

The proof is complete. □ 

Lemma 4.2. Assume (1.6), and let r > and u G L^(R''). Then is a 

well-defined function in L^iW'-). Moreover, 

(i) for a.e. x £ £^'''[m](.t) = J^^,^^ u{x + z) - u{x) dfi{z), 

(ii) ||£'^''-M||i,i(M<J) < 2||u||ii(K<i)/|^|>^dM(^;), 
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(iii) ifue L^{R'^) D BV{R'^) and r2 > n > 0, 

{{'^'"'"''[ulh^Rd) <\u\bv{R'') / \z\dfi{z) 

Jri<\z\<r2 

+ / + 2) - u||ii(Rd) d/^(^). 

J\z\>r2 



'\z\>r2 

Remark 4.3. The precise statement of (i) is the fohowing: For each a.e. repre- 
sentative of u and for a.e. x G K"^, z u{x + z) — u{x) is /^-integrable and x — >■ 
i|z|>r^(^ + ^) — u{x) d/x(z) is an a.e. representative of C^^'^[u]. This is not trivial 
since it is not immediately clear that z u{x + z) — u{x) is /^-measurable when u 
is Lebesgue measurable but not Borel measurable. 

Proof. We start by proving that 

(4.1) {z, x) e (M'^ \ {0}) X R"^ {u{x + z)- u{x)) l\z\>r is dn{z) da;-measurable 

(where we still denote by u an a.e. representative of u). It suffices to prove 
the measurability of v{z,x) :— u{x + z), since all the other terms are obviously 
measurable. Let us first assume that w is a simple function. We want to show 
that for each real interval 7, v~^{I) is an element of the (smallest) d/i(2;)da;- 

completion BQ^i\|o})xM<' of B(Rd\{o})xM<*) see the definition in Section 2.1.2. It 

suffices to prove it for u = 1^; where E G is a Lebesgue measurable set, 
i.e. when there are Bx^B^ S B^d s.t. Bi (~ E Q B2 and j^^^B^dx = 0. Now 
let Vi denote the function {z,x) — )• IbXx + z) {i = 1,2). We have to distin- 
guish between four cases according to whether / contains the values and/or 1. 
Since the proof is very similar in all cases, we do it for only one case, € / 
and 1^7. In that case v~^{I) = {{z, x) s.t. z ^0 and x + z ^ E} and v^^{I) = 
{iz,x) s.t. z and x + z ^ B,} (i = 1,2). This implies that ?;^^(/) C v-\l) C 
v^^{I), where v~^{I) € 'B(]K<i\.[o})xR<i (* = 1)2). Moreover, by the monotone con- 
vergence theorem and Fubini, 

// d/x(z) da; = sup / / IsaV-Bi {x + z) dn{z) dx 

JJ v-\l)\v-\l) n>ljR'l J\z\>l. 




sup ^ / d/i(z) / dx : n> 1> = 0. 

'b2\Bi J 

This shows that ?;~"'^(7) G '8(]K'i\{o})xR'^ hence v is measurable when u is a 
simple function. Measurability in the general case then follows since any Lebesgue 
measurable function u is the pointwise limit of simple functions. 
Now simple computations show that 



/ / \u{x + z) — u{x)\ dfj,{z) dx 

Jr'' J\z\>r 

= \\u{- + Z) -u\\Li(R<l)dll{z) <2 / d/i(z) ||M||z,i(R<i)- 

^|z|>r Jiz\>r 



<+oo by (1.6) 

By Fubini's theorem, the function z — > {u{x + z) — u(x)) l|2|>r is /i-intcgrable for 
a.e. X G R'^, and x C^'^[u]{x) := i|2|>,-"(^ + ^) ^ u(x)dfi{z) is well-defined 
in L^{R'^). This completes the proofs of (i) and (ii). When u is in addition BV, 
\\u{- + z) — ■u||Li(R£i) < |M|By(R'i) l-^l - see e.g. Lemma A. 2, and (iii) follows easily 
from the computations above. □ 

In the next result, we establish a Kato type inequality for jC^''^[A{u)]. 
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Lemma 4.4. Assume (1.4) and (1.6). Then for all u G L'^{W^), G R, r > 0, 
and6& C?°( 



sgn{u-k)C''-''[A{u)](pdx < / \A{u) ~ A{k)\ C' '''[(pjdx. 

Proof. Notice first that A{u) is by (1.4), and hence £^'''[^(u)] is well-defined 
in by Lemma 4.2. Easy computations then reveal that 

/ sgn{u- k)£''''^[A{u)](l)dx, 

sgn {u{x) — k) ^^(^(a; + z)) — ^(^(a;))^ (^(a;) diJ,{z) dx 



by Lemma 4.2, 

= / / sgn (u(a:) - k) 

JR'' J\z\>r 

{ {A{u{x + z))- A{k)) - {A{u{x)) - A(fc)) } (Pix) dij.{z) dx, 
<[ [ (\A{u{x + zj) - A{k)\ - \A{u{x)) - A{k)\) (j>{x) dij{z) dx by (2.E 

Jr'' J\z\>r ^ ' 

= 11 \A{u{x^z))- A{k)\4){x)dii{z)dx 

JRA J\z\>r 



J\z\>r 



-.1 

\A{u{x)) - A{k)\ (pix) dii{z) dx . 



Notice that the integrands above are d/x(z) dx-integrable by similar arguments as 
in the proof of Lemma 4.2. 

By the respective changes of variable {z,x) — >■ {—z,x + z) and {z,x) — ^ {—z,x), 
we find that 



1=1 I (j>{x + z)\A{u{x))- A{k)\d^i*{z)dx, 

Js.'' J\z\>r 

[ [ (l){x)\A{u{x))-A{k)\diJ,*{z)dx 

J\z\>r 



Here measure ji* in (2.9) appear because of the relabelling of z. This measure has 
the same properties as jj,. Hence we can conclude that 

f sgn{u-k)L^'''[A{u)](pdx<I- J = [ \A{u) - A{k)\ •''[(pjdx, 

and the proposition follows. □ 

The next lemma is a consequence of the Kato inequality, and it plays a key role 
in the doubling of variables arguments throughout this paper and in the uniqueness 
proof of [1, 19]. 

Lemma 4.5. Assume (1.4) and (1.6), and let u,v G L°°{Qt) n C{[0,T]; L^), 
< V e Li(M'' X (0, r)2), and r > 0. Then 

// sgn{u{y,s) - v{x,t)) 

■ (^C'^''-[A{u{;smy) - C^'^-iAivi; mix)) i,{x - y, t, s) dw < 0. 
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Proof. Note that 

sgn {u{y, s) - v{x, t)) (^A{u{y + z, s)) - A{u{y, s)fj 

- sgn {u{y, s) - v{x, t)) (^A{v{x + z, t)) - A{v{x, t)j^ 
= sgn {u{y,s) - v{x,t)) 

■ { (A{u{y + z, s)) - A{v{x + z, t))) - (^A{u{y, s)) - A{v{x, i))) } 
< \A{u{y + z, s)) - A{v{x + z, t))\ - \A{u{y, s)) - A{v{x, t))\ 

where these functions are both defined. By Lemma 4.2 (i) and an integration w.r.t. 
l|^l>rd/x(-2), we find that for all {t,s) e (0,T)2 and a.e. {x,y) e M^'', 

sgn {u - v) (r'^''- [A{u{; sMy) - C>^'^ [A{vi; tMx)) 
< [ \A{u{y + z,s))-A{v{x + z,t))\-\Mu{y,s))-A{v{x,t))\dti{z). 

J\z\>r 

After another integration, this time w.r.t. ip{x — y,t, s) dw, we then get that 
sgn{u-v) [jC'''^[A{u{;smy) - jC^'^lAivi-Mi^)) ^dw 
< / \A{u{y + z,s))-A{v{x + z,t))\ip{x-y,t,s)diJ,{z)dw 

JjQl J\z\>r 

- II I \A{u{y,s))- A{v{x,t))\il{x-y,t,s)d^l{z)dw, 

JjQl J\z\>r 

=:I + J. 

Note that the dfj.{z) dw-measurability of the integrands above follow from argu- 
ments like the ones used in the proof of (4.1). Integrability of I and J then follows 
since ||A(u)||c([o,t];L1) < ll'f'||c([o,'i'];Li) {A is Lipschitz-continuous and at 

0) and by Fubini (note the convolution integrals in x and y) 

I,J< (m(w)||c([0,T];I,i) + m(w)||c([0,T];ii))||V'llLi(K'^x(0,T)2) / dn{x). 
^ ^ J\z\>r 

To conclude, we change variables, (z, x, t, y, s) — )• {—z, x + z,t,y + z, s) in I and 
[z, X, t, y, s) — )■ (— z, X, t, y, s) in J, to find that 

1= 11 I \A{u{y, s)) - A{v{x, t))\ 'il){x + z-{y + z), t, s) dn*{z) dw, 

JJQt J\z\>r 

J = - / \A{u{y,s))-A{v{x,t))\ilj{x-y,t,s)diJ,*{z)dw. 

JJQt J\z\>r 

It follows that 7 + J = and the proof is complete. □ 
Lemma 4.6. Under the assumptions of Lemma 3.1, 

1= [[ \A{v{x,t))-A{u{y,s))\Cf[^'''{x,t,;s)]{y)dw<C, [ \z\^dfx{z), 

JJQ% JO<\z\<r 

where > does not depend on r > 0. 
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Proof. Easy computations show that 

= 9s{t -s) e^{x -y- z)-e^{x-y) + z- DOe{x - y) l|^|<i dfi*{z) 

Jo<\z\<r 

= es{t -s) j e,{x -y + z)-e^{x-y)-z- De,{x - y) l|,|<i Aii{z) 

= es{t- s)me,]{x-y), 

and by Fubini (there are again convolution integrals in /!), 

I< [[ \A{u{y,s))-A{v{x,t))\es{t-s)\jafMx-y)\dw 
JJQl 

< (\\A{u)\\LHQ^) + \\A{v)\\LHQ,))\\esC^AOe]\\LHR<i^^) 

< T\\A'\\L^(\\u\\ci[0,T];L^) + ||^;||c([0,T];Li)) II [^e] ||li(R<^+i) • 

By classical properties of approximate units (Lemma A.l in the appendix) and 
Lemma 4.1, 

=1 

< ^ll^e||w2,i(Rrf) / \z\'^d^{z), 
^ Jo<\z\<r 

and the proof is complete since 9^ e C^{R'^) in (3.1) does not depend on r > 0. 

□ 

5. Proofs of the main results 

The proofs of this section use the so-called doubling of variables technique of 
Kruzhkov [45] (see also [1, 19] for nonlocal equations) along with ideas from [46]. 
It consists in considering u as a function of the new variables {y, s) and using the 
approximate units (p'^'^ in (3.1) as test functions. For brevity, we do not specify the 
variables of u = u{y,s), v = v{x,t) and (j)'^'^ = (j)'^'^ {x,t,y, s) when the context is 
clear. Moreover, the Lebesgue measure da; dt dy ds is denoted by dw. 

5.1. Proof of Lemma 3.1. Let {x,t) e Qt be fixed and u = u{y,s), k = v{x,t), 
and (piy.s) := (jf'^{x,t,y,s). The entropy inequality for u (see (2.10)) then takes 
the form 

/ \u-v\ dscl>''^ + (qf{u,v) + \A{u) - A{v)\ bf) ■ Dy^^'^dyds 

Jqt ^ ' 

+ / \A{u)- A{y)\L^'W^\x,t,-,s)\{y)dyds 

Jqt 

+ [ sgn{u- v) C^''^[A{u{;smy)cl>''^ dyds 
Jqt 

- [ \u{y, T) - v{x, t)\ (t>''\x, t, y, T) dy 

+ / \uo{y)-v{x,t)\(l)'^\x,t,y,Q)dy>Q. 
>/m<* 

We integrate this inequality w.r.t. {x,t) G Qt, noting that Qf in (2.6) is symmetric, 
and that dscj)^'^ = -dscj)"'^ and Dy(j)'''^ = by (3.1). Consequently we find 
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that 



lie 



(5.1) 



\A{u)-A{v)\£f[<P''\x,t,-,s)]{y)dw 



+ [[ sgn{u-v)C^'^[A{u{;s))]iy) dw 
JjQl 

[f \u{y,T) - v{x,t)\ct>''\x,t,y,T)dxdtdy 
■ ff \uo{y)-v{x,t)\cl)''\x,t,y,0)dxdtdy>O. 



'QtxI 

Note that the terms in the inequahty above are well-defined since they are all 
essentially of the form of convolution integrals of i^-functions. See Lemmas 4.1 
and 4.2 and the discussions in the proofs of Lemmas 4.5 and 4.6 for more details. 
A classical computation from [46] reveals that 

h + h- II \u{y,s)-v{x,T)\(t>'^\x,T,y,s)dxdyds 



ff 

JJis.<'xQt 



(5.2) 



^ rr 



■ // \u{y,s)-vo{x)\(j)'''\x,0,y,s)dxdyds 

JJm.''xQt 



< -\\u{T) - u(r)||ii(Rd) + ||uo - 'folUi(R<i) 
+ eC^ |Mo|By(Rd) + 2w„((5) V uiv{5), 

where Cg is as in Lemma 3.1. For the readers convenience we sketch the proof 
in Section B.l in the appendix. Lemma 3.1 now follows from (5.1) and the above 

estimates on li and I5. 

5.2. Proof of Theorem 3.3. We adapt ideas from [36] to the current entropy 
solution setting by considering the region where A' > B' and its complementary. 
Let E± be sets satisfying 



(5.3) 

For all w G R, we define 



E± e Bk; 

U±£'=^ = K and r]±E±=[ 
_R\supp(^'-B')'^ C E±. 



A±{u) := r^'(r)lB±(r)dT, 

^0 

^^■^^ B±{u) := r B'{T)lEjT)dT, 

Jo 

C±{u) ■.= ±{A±{u)-B±{u)). 
We will need the following two technical lemmas. 

Lemma 5.1. Under the assumptions of Theorem 3.3, 

(i) A = A++A_ andB = B+ + S_; 

(ii) A±,B±,C± satisfy (1.4); 
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(iii) J2± |C±(w)|l1(0,T;BV) < 11^' - -B'||loo(r) \u\l^o,T;BV); 

(iv) for all ^ e M'' \ {0}, 

^ \\C±{u{- + z, •)) - C±H|Ui(Q,) < \\A' - i?'|Uoc(M) \\u{- + z,-)- wIUhq.)- 

± 

We prove these Lemmas after the proof of Theorem 3.3. 

Proof of Theorem 3.3. Let us divide the proof into several steps. 

1. We argue as in the beginning of the proof of Lemma 3.1 changing the roles of u 
and V. We fix (y, s) and take k = u{y,s) and (f)^'^ = <p'^'^{x,t,y,s) in the entropy 
inequality for v = v{x, t) to find that 

jj^^ \v - u\ dt^''^ + [qg{v, u) + \B{v) - B{u)\ bf) ■ D^^''^ dw 

+ [[ \B{v)-B{u)\Cf[cl>'''{;t,y,s)]{x)dw 

+ ff sgn{v-u)C'''''[B{v{-,t))\{x)(l>'^^ dw 
JJqI 

- ff \v{x,T)-u{y,s)\(l)'^\x,T,y,s)dxdyds 

JJM<^xQt 

+ ff \vo{x)-u{y,s)\^''^{x,0,y,s)dxdyds>0. 

JJr''xQt 

Then we add this inequality and inequality (3.3) in Lemma 3.1, 
\\u{T)-v{T)\\L^^<i^ 

< \\uo- vqWli (jg^d) +eCg\uo\BV{R'') +2u>u{S)V u)y{6) 

+ // {qg - qf){v,u) ■ D^(f)''^ dw 
JJQl 



(5.5) 



+ ff \B{v) - B{u)\ Cf [#-^(-. y, t. s)]{x) dw 

V ' 

=:/2 

+ ff \A{v))-A{u)\Cf {x,t,;s)]{y)dw 



+ 11^^ {\B{v)-B{u)\-\A{v)-A{u)\)bf -D^f^'dw 



ff ^ sgn {v - u) (^C^^'^iBivi-Mix) - C'^'^^iAiui; s))](y)) 0^-^ dw, 



where r,e > 0, < S < T, and Cg > only depends on the kernel 64 from (3.2). 
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2. It is standard to estimate Ii (cf. e.g. [25, 46] ), and I2 + 12 can be estimated by 
Lemma 4.6, 



(5.6) 
(5.7) 



h < IwoIbv(r^) Tess-supu |/' - g'\, 
h + I'2<C, [ \z\^dn{z), 

Jo<\z\<r 



where does not depend on r > 0. For the sake of completeness, the proof of 
(5.6) is given in Appendix B Section B.3. Now we focus on and I4. 

3. Cutting w.r.t. E±. We spUt and I a into four new terms using the sets E±, 
see (5.3)-(5.4). By Lemma 5.1 (i), Ia can be written as 

By Lemma 5.1 (ii), we can apply twice Lemma 4.5 with and A- instead of A, 
followed by the definitions of C±, see (5.4), to show that 



h < 



{v-u)C^^^ B+{u{;s))-A+{u{;s)) {y) 4>^^' Aw 



(5.8) 



// sgn 
JJQI 

+ ff sgn{v-u)C>''''\B_{v{-,t))-A_{v{-,t))\{x)(t)''^dw 

JJqI l j 

= 11^^ sgn {u - v) Z:"'^ [C+{u{; s))] (y) ^''^ dw 
+ ff sgu{v-u)£^'''^[C-{v{-,t))]{x)<j)''^ dw 



Note that it is crucial to have u in the first term and v in the second - otherwise 
we will not be able to apply the Kato inequality later on! 

We now consider 73. By (2.8), Lemma 5.1 (i)-(ii), the formula Dxcjf'^ = —DyCj)'^'^, 
and the definitions 1?+ = Dy and £>_ = D^, it follows that 



(\B{v)-B{u)\-\A{v)-A{u)\)d.. 

= sgn (m - w) { {A{u) - B{u)) - {A{v) - B{v)) } Dycj^''^ 

= ^sgn (u - t;) { ± {A±{u) - B±{u)) T {A±{v) - B±{v)) } D±4>' 

= J2\C±{u)-C±iv)\D±<l>^''. 



We can then rewrite I3 as 
(5.9) 



± JJqt 



4. Cutting w.r.t. z. We decompose jC'^''^ into two new terms using a new cutting 
parameter n > r. Let = /Ui + ^^\^^^^^^ for 
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and note that by (2.5), £'^''~ = C^'^^'^ + C^^'^'K Then 

1+= [[ sgn(w-t;)£'^i''-[C+K-,s))](y)<^^'Mw; 
JJQI 



(5.10) 



[[ sgn {u - v) L^^"-^ [C+(u(-, s))\{y) 0^'^ Aw. 
JJQ% 



Since C+ satisfies (1.4) by Lemma 5.1 (ii) and /ii clearly satisfies (1.6), we can 
apply the Kato type inequality in Lemma 4.4 (with k = v{x,t) and A = C+) to 
show that 

I+=[ [ sgn{u{y,s)-v{x,t))C^'''^[C+{u{;smy)<P''^dydsdxdt 
Jqt Jqt 

< f f \C+{u{y,s)) - C+{v{x,t))\C''''''[cj>''\x,t,;s)]{y)dydsdxdt. 
jQt Jqt 

Adding in the form (5.9) then gives 
(5.11) 

4+4<ll^^ \C+{u) - C+{v)\ (bf ■ Dy^^^' + C'^'^^^[cP^'\x,t,;s)]{y)) dw. 

Now easy computations show that 
Dy<f>''' = -es{t-s)De,ix-y), JC^^^^'''[<f>'''{x,t,■,s)]{y)=es{t-s)JC>'^'^[e,]{x-y). 
Hence by adding and subtracting z ■ D9t{x — y), we get that 

= 9s{t -s) [ 0,{x -y + z)-9e{x-y)-z- D0,{x - y) dij{z) 

Jr<\z\<ri 

+ 9s{t - s) D9,{x - y) ■ -6(f + / zd^l{z)] , 

V Jr<\z\<ri J 



=sgn (n -1) A(ivr)< 1^1 <ri VI ^ <i;i(^) 



where the last equality comes from (2.4) and the change of variable z —z. We 
insert (5.12) into (5.11) and combine the resulting inequality with (5.10), 

It + It< 

[[ \C+{u)-C+{v)\ 
JJQI 

■ 0s{t -s) 9^{x -y + z)- 9^(x - y) - z ■ D9^{x - y) d/x(^;) dw 

Jr<\z\<ri 

(5.13) + jj^^ \C+{u) - C+{v)\ 

■ 9s{t — s) D9^{x — y) ■ sgn {n — 1) / zdfj,{z)dw 

J ri A(lVr)<|0|<riVl 

+ [[ sgn(u-z;)£^''-i[C+(«(-,s))](y)</>^'^d«; 
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Similar arguments show that we can bound /g + (see (5.8)-(5.9) ) as foUows, 

l3+l4< 

[[ \C.{v)-C.{u)\ 
JJol 



■ Os(t -s) e^ix - y - z) - Oeix - y) + z ■ D9^{x - y) dAt(z) dw 

Jr<\z\<ri 

\C.{v) - C.{u)\ 
(5.14) -^-^Qi- 

■ 6s{t — s) D6^{x — y) ■ sgn(ri — 1) / zAiJ,{z) 6.w 



sgn {v - u) £"''■1 [C-{v{-,t))]{x) (j)^^^ Aw 



<JJq2 sgn{v-u)Cf-^i[C-{u{-,s))]{y)4><''^dw by Lemma 4.5 

=: Ji + J2 "I" "^3 • 



5. n BV -regularity. It remains to estimate for i = 1. . . . , 3 in (5.13) (5.14). 
For and , we use Fubini and integrate by parts to take advantage of the BV- 
regularity of the entropy solution u. After some computations given in Appendix B 
(see Lemma B.l), we find that 

\jt\<^ I \DeMx [ kPd/z(z)|C±(«)Ul(o,T;By), 
JW^ Jr<\z\<r, 



\4\< 



j 

7riA(lVr)< 



^d/x(z) 



riA(lVr)<|z|<riVl 



|C'±(w)|Li(o,T;By), 



and hence 



V(Ji± + 4) <^ [ Wdx [ \Z\' dl,{z) V |C±(w)Ui(o,T;W) 

_l_ jRd Jr<\z\<ri _|_ 



+ 



/ 



zdiJ,{z) 



\C±{u)\li{0,T;BV)- 



' riA{lVr)<\z\<riVl 

By Lemma 5.1 (iii) and a priori estimates for u, cf. (2.11), we see that 
V(J± + J^) < 1 / \DeM^ \u\lho,T;BV) [ \zfdt,{z) \\A' - B'^o 

_l_ Jm-i ^ V 'Jr<\z\<ri 

<l«olBV(Kd)r 



(5.15) 



Ml^(0,T;BV) 



^I"oIbv(r<') ^ 



riA(lVr)<|z|<riVl 



;d/z(2;) 



Let us now estimate Jg" in (5.13). Easy computations (see the proofs of Lem- 
mas 4.5-4.6 ) show that 



4 < \\Os0e\\LHM^+^) ||£'^''-^[C+(n)]|Ui(Q,). 
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By Lemma A.l, \\0s 9e\\L^(Rd+i) = ||^5||z,i(r) ||6'£||Li(Rd) = 1, and then we can use 
Lemma 4.2 (iii) to find that 



4< 



f if \C+{u{-,s))\BviR'^)\z\dfi{z) 

Jo \^ri<|2|<r2 

+ / \\C+{u{- + z,s))-C\{u{-,s))\\mR<i)dn{z)]ds 

J\z\>r2 J 

for all ra > n. Since C+(tt) G n C{[0,T];L^) n L°°{0,T; BV), {z,s) 
||C+('u(- + z, ,s)) — C+(u(-, s))||^i(][{d-) and s ^ |C+(m(-, s))|5y(][{d) are continuous and 
lower semi-continuous functions respectively and hence Borel measurable. They are 
thus diJ.{z) ds-measurable, and we may change the order of the integration to find 



4< 



/ \z\dlJ,{z)\C+{u)\Li(^o,T;BV) 
Jri<\z\<r2 

+ f \\c+iu{- + z, •)) - c+(n)|Ui(Q,) d^x{z). 

J\z\>ro 



'\z\>r2 

We get a similar estimates for and find by Lemma 5.1 (iii)-(iv) and (2.11) that 



< / \z\dlJ,{z)^\C±{u)\Li(0,T;BV) 
_|_ Jri<\z\<r2 ± 

+ f E ll^±("(- + ^' •)) - C±{u)hHQ^) df^iz), 

J\z\>r2 ± 

< |wo|By(R<i) T / 1^1 dn{z) \\A' - B'||ioo(K<i) 

J ri<\z\<r2 

(5.16) + / \\u{- + z, •)) - wIUhq,) dM(^) \\A' - B'U^^^.y 

J\z\>r2 V ' 

<T||«o(-+z)-«ollil(Kd) 

The last inequality (under the bracket) comes from (2.12) applied to the solu- 
tion u{- + z, •) of (1.1) with initial data uo{- + z). 

6. Conclusion. By (5.8)-(5.9) and (5.13)-(5.14), h+h < E± Ei=i ^t- Therefore 
we may estimate (5.5) by (5.6)-(5.7) and (5.15)-(5.16). For all r2 > ri > r > 0, 
e > 0, and T > S > 0, we find that 

\\u{T) - w(T)||ii(R<i) 

< \\uo - vo||z,i(iR<i) + |uo|BV(i:<i) ^'ess-supR |/' - g'\ 



(5.17) 



+ eCg\uo\BV{R^)+'^^u{S)\/ujy{6) + Ce / \zfdn{z) 

Jo<\z\<r 

+ 1 / \D0d\dx\uo\Bvim^)T f \z\^ dfi{z)\\A' - B'U 

JR'' Jr<\z\<ri 



\\A'-B'\\^. 



/ zdii{z) 

J riA{lVr)<\z\<riVl 

+ \uo\bv{r<^) T / 1^1 dii{z) \\A' - B'lli^fRd) 

J ri<\z\<r2 

+ T \\uo{- +z)- uoWmQT) dn{z) \\A' - B'||i«>(R.), 

J\z\>r2 



'\z\>r2 

where Ce > does not depend on r > 0. 
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To finish, we first pass to tlie limit as r — >■ in (5.17). By tlie dominated 
convergence theorem, the result is equivalent to setting r = in each term, and in 
particular the term Ce Jo^^z\<r l-^P'^A*!-^) vanishes. Secondly, we pass to the limit 
as 5 — to get rid of the term 2lVu{6) V uiy{S). Finally, wc optimize the remaining 
terms w.r.t. e > by using the formula min£>o (ea + ^) = 2\/ab (for a,b > 0). 
This gives us the following continuous dependence estimate: F or all r2 > ri > 0, 

v\\c([o,T];L^) < \\uo - vo\\L^(s,d) + \uo\BV{s.d) T ess-swp^ \g' - /'I 
+ 2JJq / \D9d\dx\uo\ly^^,.T [ |z|2dM^)||A'-iJ'|U=o(R) 

y ^ -/Md ^ ' Jo<\z\<ri 

+ \uo\BViR-i)T / \z\dn{z)\\A' - B'Wlo 

J n <lzl<r9 

+ |woUy(iR<*) ^ 
+ T I \\uo{- + z)- Uo|Ui(R<i) dfiiz) \\A' - B'||ic»(R), 

J\z\>r2 

where 64 is an arbitrary approximate unit (3.2) and Cg = 2 J^^ \x\ Od{x) dx by 
Lemma 3.1. 

Let 0d = On where {On}neN is a sequence of kernels s.t. 0„ satisfies (3.2), 0„ — >• 
Wd~^l|.|<i in and f^^ \D6n\ dx Wd~^|l|.|<i|si/(Rd). Here uid is the volume of 
the unit ball in M'*. Note that the i?T/^-semi-norm of the indicator function of the 
unit ball is equaled to the surface area of the unit sphere, i.e. |1|.|<i|bv'(K'') = <^<^d- 
Moreover, we have 



(5.18) 



/0<|z|<ri 

= (R) 

ri<\z\<r2 

zd^I{z) WA'-B'W 

riAl<\z\<riVl 



f |a;||^„(a;)| da; — ^ — f \x\dx = 

JM-i J\x\<l 



d 



J\x\<l d+1 

The proof of (3.4) is then complete after passing to the limit as n ^ +00 in 
(5.18). □ 

Let us now prove Lemma 5.1. 

Proof of Lemma 5.1. The proofs of (i)-(ii) are easy and left to the reader. Let us 
prove (iii)-(iv). We differentiating (5.4), = ±{A'^-B'^) = ±{A'-B')1e± a.e., 
and use (5.3) to see that 



(5.19) C'^ = {A' - B' 



a.e.. 



Since u{t) e L^{W^) n BV{W^) for fixed t G (0,T), there is {0„}„6n C C°°(R'^) n 
W^'^(R''^) s.t. (pn u{t) in and J^^ |D0„|dx — > \u(t)\Bv see e.g. Theorem 
2 in Section 5.2 of [33]. Since C± satisfy (1.4), it follows that C±((/)„) C±{u{t)) 
in L^. Moreover, by lower semi-continuity of the BT^-semi-norm, 

V \C±{u{t))\Bvm ^ El^i^f [ \DC±{M\ dx < liminf / V \DC±{<f>n)\ dx. 

Since <^ri is smooth and C± is Lipschitz-continuous, we can use the chain rule 
and (5.19) to show that 

53 |DC±(<^„)| = |A' - B'|((/.„) < P' - B'IUoo(R) \Dc^n\ a.e.. 
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In the limit as n — > +00, one thus gets 

^\C±{u{t))\BV(R'i) 
± 

< P'-S'llioo(K) hm / |£>0„|da;= p'-B'||i<»(R)|w(f)|Bv(Ma). 

An integration w.r.t. t £ (0,T) now completes the proof of (iii). 

We now prove (iv). Since C± is Lipschitz-continuous, it can be written as the 
integral of its derivative. This implies for a.e. {x, t) G Qt, 

Y,\C±{u{x + z,t))-C±{u{x,t))\ 
± 

= \u{x + z,t) - u{x,t)\y2 \ ( C'±{{'^-T)u{x,t) + Tu{x + z,t))dT 

<\u{x + z,t)~u{x,t)\ \A'-B'\{{1 — t) u{x, t)+Tu{x + z,t))AT by (5.19), 
< - \u{x + z.i)-u{x,t)\. 

The proof of (iv) is complete by integrating w.r.t. [x, t) € Qt- □ 

5.3. Proof of Theorem 3.4. We argue step by step as in the proof of Theorem 3.3. 

This time, E± are taken such as 



(5.20) 



(5.21) 



E± G Sr; 

U±£;± = M and n± E± = 0; 
^M'*\supp(/x- !/)'F C E±. 

Let fjb± and v± denote the restrictions of and u to E±. It is clear that 

,M±ji^±5 and ± {n± — u±) all satisfy (1.6). 
Proof of Theorem 3.4- 

1. We apply Lemma 3.1 with A = B, but different Levy measures /j. and v, along 
with the entropy inequality for v to show that for all r, e > 0, < 5 < T 

||u(T)-u(T)||ii(R.) 

< ||wo - ^^o||z,i(K<i) + ^C§\uo\bv(R'^) +2uJu{6)V u>y{6) 
+ // iqg-qf){v,u)-D^(l)''^dw 



(5.22) 



If \A{v)- A{u) I £f [cl>''\; y, t, s)] (x) dw 
ff \A{v)) - A{u)\Cf [,^^'\x, t,;s)]{y) dw 
11^^ \A{v)-A{u)\ [bf -bf) ■D,ct>'''dw 



+ 



ff ^ sgn {v - u) [A{v{;tmx) - /:''''^[A(i.(., sMy)) 4,^'' dw, 
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where > does not depend on r > 0. Except for I3 and I4, the other terms 
were estimated in the proof of Theorem 3.3. 

2. Cutting w.r.t. E±. We use the notation introduced in (5.20). We apply 
Lemma 4.5 twice with u+ and /X- instead of /x, along with linearity of £^''^ in 
/i, see (2.3), to see that 

< [[ sgn{v-u)(c^^'^[A{u{-,s))]{y)-C^+'^[A{u{;s))]{y)) <t>'^' dw 

+ 11 ^ sgn(t;-n)(£^-''-[A«,i))](a;)-£''-'^[A«,i))](a;)) <^'^' dw 

= /7 sgn(u-t;)/:^+-''+''-[A(u(-,s))](y)(/)^'''dw; 

+ sgn(i; -«)£-('*— "-)'''[A(t;(-,i))](a;)<^^'^d«; 
JJQI 

(5.23) -:/4+ + /4-. 

Again, it is crucial to have u in and in in order to use Kato's inequality 
later on. 

Let us now consider /s. By (2.4) and (2.9), 6^ and /x* are linear w.r.t /U. Easy 
computations using (5.21) then leads to 

where £>+ = and £>_ = £>x) and hence 

^3 = E // - ^(^)l • ^±<^''' d«; =: /3+ + /s". 

3. Cutting w.r.t. z. The computations of this step are similar to the ones in the 
proof of Theorem 3.3. For the reader's convenience, we estimate +1^ , the terms 
that was left to the reader in the preceding proof. 

For any measure jl we let jli = /U|o<|^|<ri and write fb = p,i + ^^\^^^^^^ for ri > r. 
Then 

I4 < [[ sgn(t;-u)£-(''--^-)i''-[yl(u(-,t))](a;)0^'^dw 

V ' 

+ [[ sgn{v-u)C-''^'—'-^'''^[A{v{-,t))]{x)(l)'''^dw. 
JJQt 

Recall that —{p- — is a positive Levy measure by (5.21), so we can apply 
Lemma 4.4 with — (/i- — '^-)i instead of /i and k = u{y, s) to find that 

I^< ff \Aiv) - Aiu)\C-^^—'-^'^''-[cl>^'\;t,y,s)]{x)dw 

and 

jj^^ \A{v)-A{u)\ (b.-f''— '^-)* -D.^^'* + '^-)^''"[,^^'*(.,i,y,s)](a;))dw;. 
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Easy computations then leads to 

= es{t -s) f 9,{x -y-z)- e,{x - y) d(z/_ - ii-){z), 

J r< |2:| <ri 

and we can rewrite the nonlocal operator as follows, 
&-(^— )* ■ D,<t>'^' + -^-y-'^ [(/.^'^(•, t, y, s)] {x) 

= 9s{t -s) f e^{x -y- z)-e,{x-y) + z- DO^ix - y) d(i/_ - M-)(^) 

Jr<\z\<r-^ 

- 9s{t - s) D9,{x - y) ■ (-b;^>^—-y + j 

\ Jr< 



r<\z\<ri J 

V 

V 

=sgn (ri-l)X.^^(^y^,<|^|<^^y^zd(i/_-Ai_)(0) 

Compare this expression with (5.12) that appear when and I4 are considered. 
We add the different estimates and find that for all n > r, 



< ff \A{u)-A{v)\9s{t-s) 

JjQl 

■ / 9e{x -y-z)- 9e{x -y) + z- D9^{x - y) d{v- - H-){z) Aw 

J r< \z \ <ri 

- // \A{u) - A{v)\ 9s{t - s) D9e{x - y) 

JjQl 

•sgn(ri — 1) / zd{i'- — ii-){z)dw 

J riA{l\/r)<\z\<ri\/l 

+ ff sgn{v-u)JC-'^^'—'-y'''[A{v{■,t))]{x)(l)''^dw 
JjQl 



<ffQ2 sgn{v-u) C "-'■'■i [A(tt(-,s))](2/) du) by Lemma 4.5 

— Jl + ^^2 + "^3 • 

Similar arguments also lead to 

4 + It 

< ff \A{u)-A{v)\9s{t-s) 

■ I 9e{x -y + z)- 9e{x -y)- z- D9e{x - y) d(/u+ - iy+){z) dw 

J r<\z\<ri 

+ ff \A{u) - A{v)\ 9s{t - s) DUx - y) 

■ sgn in - 1) / z d(/i+ - iy+){z) dw 



'riA(lVT-)<|2|<riVl 

+ ff Sgn {u - v) [A{u{-, 0^'^ dw, 

=: J+ + ... + J3+. 
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4. n BV -regularity. We estimate {i = 1, . . . , 3). By (B.l) of Lemma B.l and 
(2.11), it follows that 

E-^i^^^ / \DeMx\uo\Bv(M'')T\\A'h^^^) [ |^|2dV±(/x±-i/±)(^). 

_|_ JR"^ Jr<\z\<ri _|_ 



by (5.21) 



Note now that X^_|_(m± — ^^±) = M ~ hence 

VJ2±= //' |A(w)-A(t;)|(?5(t-6-)i?^,(x-y) 

•sgn(ri — 1) / zd{fi — i'){z)dw. 

J riA{lVr)<\z\<riVl 

An other apphcation of (B.2) of Lemma B.l and (2.11), can be used to see that 
^ J2=^ < |uo|sy(Md)T||A'||ioo(R) / zd{i^- iy){z) . 

_|_ ^riA(lVr)<|z|<riVl 

Finally, we use again Lemma 4.2 (iii) and (2.12) to show that 

< \uo\BV(R'i)T\\A'\\i^oo(Rci) / \z\d\ii- fliz) 

j- Jri<\z\<r2 

+ rp'||ioo(R) / ||Mo(- + ^)-'«o||Li(R^)d|/X-I/|(z). 

J\z\>r2 

5. Conclusion. The rest of the proof is the same as for Theorem 3.3; i.e. we use 
the estimates on jf to estimate I3 + 14 < J2i^i J2± (5-22) and pass to hmit 

and/or optimizes w.r.t. the parameters r, e, 5 > 0. The proof is complete. □ 

Appendix A. Properties of approximate units and sy-PUNCTioNS 
The results in this section are quite classical, see e.g. [13, 33, 53]. 

Lemma A.l. Assume (1.6) and let (j)^'^ be defined as in Lemma 3.1. Then for 

all e, (5 > 

(i) II^(S||li(E) = ll^'cllLi(R'i) = 1, 

(ii) ^^,\De,\dx<\\\Dh\\mm- 

Lemma A.2. Let u e L'^{W^) n BV{W^). Then for all h G W^, 



/ \u{x-\-h)-u{x)\dx<\h\\u\Bv{R'i)- 



Lemma A.3. Let u e L^{W^) n BF(M<^) and g : R ^- M*^ he Lipschitz-continuous. 
Then for all </> G C° 



(A.1) 

Moreover, if 77 : 



q{u) ■ D(j)dx 



< ess-supg |g'| / \(j>\d\Du\{x). 



is Lipschitz-continuous, then r]{u) G 
\Dv{n)\ < \\v'\\L^m\Du\. 



with 



Let us give a short proof of the last lemma for the reader's convenience. Re- 
member that I • I denotes the total variation measure of a Radon measure (see the 
definition in (2.1)). 
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Proof. Let us prove (A.l). For all n G N, define := u*0^ G C°°{R'^)nW^''^{M.'^) 
where Oi. is the approximate unit in (3.1). Also assume 9^ in (3.2) is even. By 
classical results on approximate units, (f)n ^ u in as n ^ +oo. Moreover, for 
all X e W^, 



\D<Pn{x)\ 



9xix-y)dDu{y) 
<f 0.{x-y)d\Du{y)\= f Oi{y-x)d\Du{y)\, 

J'R<' " JlRd " 



since 9i is nonnegative and even. It then follows that 



In- 



g(0„) • D(j)dx 

4>q'{(l)n) ■ D(j)n dx 
< ess-supjj l^'l / \D<p., 



< ess-supg \q' 



''I / 1-^(^)1 / 



by integration by parts and the chain rule, 

da;, 

Oi {y-x) d\Du\{y)dx. 



Since \Du\ is cr-finite as finite Radon measure, the product measure dxd\Du\{y) is 
well-defined and Fubini's theorem applies. After changing the order of the integra- 
tion, we find that 

In < ess-supR|g'| / \(f)\ * 62^{y)d\Du\{y). 

Again by classical results on approximate units, * 9i ^ \<p\ uniformly on R'^, 
and hence 

limsup/„ < ess-supgl^'l / |0| d|£>u|(y). 

Finally since q{(l>n) converges to q{u) in Lj^^, we also have that 

lim /„ = / q{u) ■ D(j)dx , 

and the proof of (A.l) is complete. 

To prove the assertion for r]{u), we repeat the above arguments with r] instead 
of q. The result is that for all (/) e C;?°(R'^,M'^), 

/ r?(w)div</)dx < ||r?'|Uoo(K) /" \(l>\d\Du\{y), 

and we may conclude from Riesz's representation theorem [33, Section 1.8 Theo- 
rem 1] that r]{u) e BViW^). Moreover, by (2.1) , we see that for all B G B^^d 

\DT]iu)\iB) < Wv'h^m mi{\Du\{U) s.t. U open, B C U} . 

Finally, by classical regularity properties of Radon (Borel) measures [33, 53], the 
infimum above is equaled to \Du\{B) and the proof is complete. □ 

Appendix B. Technical results and computations 
B.l. Proof of (5.2). We start by estimating I4. Since 

-\u{y,T)-v{x,t)\ < \u{y,T) - u{x,T)\ - \u{x,T) - v{x,T)\ + \v{x,T) - v{x,t)\, 
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and (j)'^'^ is nonnegative, 



h = - II \u{y,T) - v{x,t)\ct>'-'{x,t,y,T)dxdtdy 

< 



// 

JJQt> 

II \u{x,T)-v{x,T)\(j)'^\x,t,y,T)dxdtdy 
+ f[ \u{y,T) - u{x,T)\(f>''\x,t,y,T)dxdtdy 

JJQtX'R'' 

+ // \v{x,T) -v{x,t)\(l)''-^^{x,t,y,T)dxdtdy 
=: J1 + J2 + J3. 

By Lemma A.l, 

Ji = - [ \u{x,T) - v{x,T)\ [ 0s{t-T) [ e,{x-y)dydtdx 
JR'' Jo Jm.'' 

^ V ' 

=1 

= -h(r)-^;(r)|Ui(R.) [h{T)dT. 

By the change of variables {x, t, y) {x, t,x — y), 

J2= I es{t-T)ee{y) [ \u{x-y,T)-u{x,T)\dxdydt, 

JQt •'r<* 

(r)lBy(R.) [ es{t-T)\y\e,{y)dydt by Lemma A.2, 

Jqt 

f h{T)dT [ \y\^~9Jy)dy by(2.11), 
J -3^ jRd Ve/ 

\y\Odiy)dy luolsviR-') /^^iWd-r. 



< u 



< 



Finally, we see that 



J3= [ \v{x,t)-v{x,T)\Os{t-T) [ 9,{x-y)dydxdt, 
<[ es{t-T) [ \v{x,t) -v{x,T)\dxdt since supp C [-5,(5] by (3.2), 

,{d) j\ h{T)dT. 



nO 



We conclude that 



Starting from 

-\u{y, s) - v{x,T) \ < \u{y,s) - u{x, s)\ - \u{x,T) - v{x,T) \ + \u{x, s) - u{x,T)\, 
similar reasoning leads to 
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Hence, h + I^ < -||it(T)-'u(r)||ii(Rd) + ^ e |Mo|By(Rd)+w„(<5)Vw„(5) when ^ > 1 
and we finish the proof of (5.2) by estimating + with similar arguments. 

B.2. Estimates of the singular nonlocal parts. 

Lemma B.l. Assume (1.4) and (1.6). Let u,v G L°°{Qt) n C{[0,T]; L^) n 
L°°{0,T;BV), (jf^^ he as in Lemma 3.1, and n > r > 0. Then 



[[ \A{v{x,t))-Aiuiy,s))\ 
JJqi 

6e{x — y± z) — 6f{x — y)^ z ■ D9^(x — y) diJL{z) dw 



(B.l) • es{t-s) 



r<\z\<r\ 



1 



2e 



<- f \D~ed\dx [ 



r<\z\<ri 



\zfdlJ,{z) |^(u)Ui(0,T;W)- 



and 



(B.2) 



// \A{v{x,t))-A{u{y,s))\ 
JJQI 

es{t - s) D6,{x - y) ■ sgn{rx - I) f 

JriA(l 

I 



< 



/riA(lVr)<|z|<riVl 



Vr)<|z|<riVl 
Zdll{z) \A{u)\l1(0,T;BV)- 



z diJ,{z) dw 



Proof. We start by proving (B.l) in the + case. Similar arguments give the proof 
also in the — case. From Taylor's formula with integral remainder, 

deix - y + z) - Oeix -y)- z- DOeix -y) = [ {1-t) D^O^ix -y + Tz)z- zdr. 

Jo 

Let / denote the integral in the left-hand side of (B.l). By Fubini's theorem, 
(B.3) /= // 

JJ (0,Ty Jr<\z\<rx JO 

, t)) — A{u{y, s)) \ D^9e{x ~ y + t z) z ■ zdydx dTd^{z) dtds. 



Jr<\z\<ri Jo 



[ [ \A{v{x,^ 



= :J 



The dr dii{z) dw-integrability of the integrands follows from similar arguments as 
in the proof of Lemma 4.2. Note that r]k{A{u{- , s))) = |fc - A{u{-,s))\ e BV with 

\D^k{A{u{-,s))\ < \DA{u{;s))\ 

for any /c € IR. by Lemma A. 3 and L-^ D BV regularity of A{u{-, s)). Integration by 
parts w.r.t. y (for fixed z, x, t, s), then leads to 



Ul 



/ f De,{x- 

■/R<i JR 



y + Tz)- zz- dDT]A(v(x,t))iA{u{-,s))){y) dx 



\D0,{x -y + Tz)\ d\DA{u{; s))\{y) dx. 
We change the order of integration (using Fubini) and use Lemma A.l to see that 

\J\<\z\''\A{u{s))\bv(R'^) [ \De,{x)\dx<\z\''\A{u{s))\BviR'')- [ \Ded\dx, 



CONTINUOUS DEPENDENCE ESTIMATES FOR INTEGRO-PDES 



31 



11(0 



and then from (B.3) that 
|/|<- / \D0d\dx 

f [ es{t^s){l-T)\z\^\A{u{s))\By(^K.-^dTdn{z)dtds. 

'(0,T)2 Jr<\z\<ri Jo 

Let us recaU that the integrand above is dTd/i(2) dt ds-measurable since s ^■ 
\u{s)\bv(R'') is lower semi-continuous. By Fubini we then integrate first w.r.t. t 
and use that 6s {t — s)dt <1 to see that 

\I\<- [ \D~ea\dx [\l-T)dT [ \z\^di,{z) [ \A{u{s))\BviR'')ds, 

£ ^M<i Jo Jr<\z\<ri Jo 

and the proof of (B.l) is complete. 

We prove (B.2) by similar arguments. Define 

(B.4) q{v,u) := \v — u\sgii{ri — 1) / zdiJ,{z), 

J riA(lVr)<\z\<riVl 

and note that it is Lipschitz-continuous. Again we denote by / the integral of the 
left-hand side of (B.2). By Fubini's theorem, 

(B.5)/=/7" 0g{t-s) [ [ D0,{x-y)-qiAiv{x,t),Aiu{y,s)))dydx dtds. 

Jj{0,TV J'SA jR'i 

= :,/ 

For fixed {x,t), q{A{v{x,t), ■) is Lipschitz-continuous and we may use (A.l) and 
integration by parts in y to see that 



I J| < ess-supR2 \qu\ / 0, 



{x-y)d\DA{u{-,s))\{y)dx, 



where ess-sup]g2 \qu\ denotes the Lipschitz constant of q w.r.t. its second variable 
u. Changing the order of integration, we find that 

J < |A(w(s))|bv(k'')CSS-supr2 \qu\, 

and hence by (B.5) and integrating first w.r.t. t, we get that 

(B.6) |7| < ess-supR2 |g„| / \A{u{s))\Bv(R'i) ds. 

Jo 

The proof of (B.2) is now complete since by (B.4), 



ess-sup]K2 \qu\ 



/ zdfM{z) 

^riA(lVr)<|z|<riVl 



□ 



B.3. Proof of (5.6). It remains to prove (5.6) from the proof of Theorem 3.3. 
Recall that (/ is the entropy solution to (1.1) and that Ji is defined in (5.5). We 
want to prove that 

I '■= // {Qg-Qf){v{x,t),u{y,s))-0s{t-s)D0^{x-y)dw 
JJQI 

< |Mo|By(M<i) Tess-supK |/' - fl^'l- 

Define 

q{v, u) := {qg - qf){v, u) = sgn {v - u) {{g{v) - g{u)) - {f{v) - /(u))} 
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(see Definition 2.6). Since /, 5 are Lipschitz-continuous, q will be so too, and / 
will satisfy (B.5) with A{u) = u and new flux q. Arguing word by word as in the 
preceding proof from (B.5) until (B.6) using also (2.11), we find that 

\I\ < ess-supR2 |g„| / \u{s)\BviR-') - l^*o|By(K'') I'ess-supRS |g„|. 
^0 

The proof is complete since ess-sup]^2 | qu \ ^ ess-supj^ | — by definition. 
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